Compute result from input 
using any simple method 
known in the art. 



Set n2 to n.div.2. 



i 



Set n3 to n-n2. 



I 



Compute the first partial 
solution x, using the first part 
recursively at precision n3. 



Compute y using the Frobenius operator. 



i 



Compute the error term e from x and y. 



Divide e by p A n3. 



i 



Compute correction factors u and v. 



Compute the second partial 
solution z, using the second 
part at precision n2. 

Return x + z * p A n3. 




Figure 1 : First part for precision n. 




Set n2 to n.div.2. 
Set n3 to n-n2. 



Compute the first partial 
solution x, using the second part 
recursively at precision n3. 



Compute y using the Frobenius operator. 
Update the error term e using x, y and u, v. 

j ~ 

Divide e by p A n3. 



Compute the second partial 
solution z, using the second part 
recursively at precision n2. 



i 

Return x + z ■ p A n3. 




Compute result from e 
using any simple method 
known in the art. 



Figure 2 : Second part for precision n, with error 
term e and correction factors u and v. 



